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For a periodically shaken optical lattice, effective time-reversal is investigated numerically. For 
interacting ultra-cold atoms, the scheme of [J. Phys. B 45, 021002 (2012)] involves a quasi- 
instantaneous change of both the shaking-amplitude and the sign of the interaction. As the wave 
function returns to its initial state with high probability, time-reversal is ideal to distinguish pure 
quantum dynamics from the dynamics described by statistical mixtures. 
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I. INTRODUCTION 

In the classical world our intuition is based on, there 
are many processes which cannot easily be reversed. For 
many cases a film run backwards will be easily identifi- 
able as such: a broken glass does not reassemble itself 
and jump on the table; a mixture of milk and tea does 
not unmix in such a way that the milk flows back into 
the milk bottle. However, the same is not always true 
for quantum mechanics. This paper investigates effective 
time-reversal [T] via periodic shaking (FigJT]) for ultra- 
cold bosonic atoms in an optical lattice [2HIJ7 




FIG. 1. (Colour online) Ultra-cold atoms in an optical lattice 
combined with periodic shaking. The full quantum dynamics 
are modelled by the time-dependent Bose-Hubbard Hamilto- 
nian |TJ; many aspects can be understood by the effective 
Hamiltonian jjjt. 

Tunnelling control via periodically shaking an opti- 
cal lattice has already been applied in several experi- 
mental setups. The topics investigated experimentally 
range from control of the superfluid-to-insulator transi- 
tion [5] (cf. [6l [7]) over frustrated classical magnetism [8] 
to photon-assisted tunnelling [5HPT] (cf. [T2HIS])- Re- 



search on periodically shaken systems also includes de- 
struction of tunnelling [ToT - H5] . generation of quantum 
superpositions [201 HI] , robust dynamical recurrences [22] 
and reflection- less defects (22] . 

Experiments with periodically shaken optical lat- 
tices [5] [TU1 [TTJ [21] can often be modelled via a time- 
periodic Bose-Hubbard Hamiltonian like 

oo 

H(t) = - j + 

l=-oo 
oo 

+ 7T 2^/ C t C t C t C l 

e=~oc 

oo 

+ ^ 2£h[i cos(ut + a)c\c e . (1) 

l=-oo 

The operators c»/cj annihilate/create a boson in well i; 
J describes the hopping of particles in the lattice and 2hfi 
is the driving amplitude. The on-site pair interaction is 
denoted by U. In the high-frequency limit (fiw U, 
hu 3> J), many aspects of experiments can be under- 
stood by an effective-Hamiltonian approach for which 
Eq. ([I]) is replaced by [5] 

oo 

H cf f = - Jeff ( d l d £+l + 4+1^) 

e=-co 

+ 2 E 44 £ e £ t > ( 2 ) 

t=-oo 

where the hopping J is multiplied by the Jq Bessel- 
function [23] 

Jeff = J Jo (£) • (3) 
For larger interactions, the model ^ has to be modified 

(cf. urns]). 

In experiments, systems modelled by time-periodic 
Hamiltonians can be switched on quasi-instantaneously. 
If - as in Eq. ([I]) - the Hamiltonian contains a term 
oc cos(ujt + a), it does make a difference if the cosine is 
at a minimum or maximum. Counter-intuitively, switch- 
ing it on at a maximum is the best thing to do if one 
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does not want to change the wave function (cf . [57] ) . 
For periodically shaken systems of ultra-cold atoms, sud- 
den changes of shaking have been used to cover topics 
as diverse as super Bloch oscillations [53] [55] [5J5], di- 
rected transport [30], effective magnetic fields [3T] and 
time-reversal for a Bose-Einstcin condensate in a double 
well [I]. 

The topic of this paper is to apply the time-reversal 
scheme of Ref. [T] to periodically shaken lattices. Sec- 
tion [n] motivates time-reversal as a method to distin- 
guish quantum dynamics from similar looking classical 
diffusion. The model used to describe the shaken lattice 
is introduced in Sec. Ill (and in the Appendix). Sec- 
tion HYI introduces the time- reversal scheme. The results 
are presented in Sec. [V] 



II. QUANTUM DYNAMICS VS. CLASSICAL 
DIFFUSION 



Figure [2] shows both the solution of the one- 
dimensional diffusion equation and of the one- 
dimensional Schrodinger equation for a free particle. For 
Fig. [2]both the classical diffusion equation in one dimen- 
sion 132 , 



d t g(x,t) = DB%Q(x,t) 



(4) 



where D > is the diffusion constant, and the 
Schrodinger equation for a free particle of mass m in one 
dimension, 



2m 



(5) 



have been solved for the initial conditions g(x, t — 0) = 
exp(— x 2 )/^ and \t'(x,t = 0) = exp(— 0.5 x 2 )/-^ using 
dimensionless variables (D = 1, m = 0.5, h = 1). Both 
the solution of the diffusion equation [35] and the solution 
of the Schrodinger equation |33j are known analytically. 
The rate at which the width of these distributions spreads 
is a valuable way to distinguish quantum dynamics from 
classical dynamics 34J. 

On the quantum level it is, in principle, possible to 
construct an initial wave function for which the width 
initially decreases by, e.g., numerically preparing a wave 
function for which the probability function reverses the 
motion shown in Fig. [2] (b). There is even an exam- 
ple where such a case was realised experimentally: by 
cleverly changing the wave function [35 , the quantum 
motion has been turned backwards in time for atom-field 
interaction in a cavity quantum electrodynamics exper- 
iment 36 . Time-reversal schemes on the level of the 
Gross-Pitaevskii equation can be found in Refs. [37J [35] . 



III. MODEL 

This section derives the effective Hamiltonian ^ for 
a single particle in a double well; the extension to an 




(c) 



<D 
"D 



0.5 ■ 




FIG. 2. Quantum dynamics versus classical dynamics in one 
dimension, (a) Two-dimensional projection of the solution 
g{x,t) of the classical diffusion equation Q. To increase the 
visibility, g ' 25 was plotted (in arbitrary units), (b) Two- 
dimensional projection of the probability density |$(a;,f)| 2 
of a free particle [cf. Eq. Again, (jv]/| 2 ) - 25 is plotted, 

(c) The initial conditions for both cases are chosen to be the 
same. 



optical lattice is given in the Appendix. 

For a single particle in a double well, the Hamilto- 
nian reduces to 



H(t) = -j(c\c 2 + clc 1 



+h/i cos(ujt + a) C 



(G) 



one thus has to solve the corresponding Schrodinger equa- 
tion: 

id t ipi (t) = -[i cos(ujt + a)ipi (t) - ^ip 2 (*) (7) 

idti>2{t) = +ficos(ujt + a)ip2(t) - T ipi(t) (8) 

n 
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The ansatz 



(t) = ipi (t) exp i—sm(tjjt + a) 



4>2{t) = <p2{t) exp — i— sin(o;< + a) 



leads to 

id t tpx(t) 



(p 2 (t)e-xp 
<^i(i)exp 



2 M • , 
% — sin(wj + a j 

2/" • / 
— i — sin(wi + a) 
to 



(9) 
(10) 

(11) 
(12) 



In order to solve this in the high-frequency limit Hu) 3> 
J (which may already be reached for hu; ss 10 J), the 
expansion in Bessel functions [35] is particularly useful 



a iz sin(u;i-t-. 



(13) 



For high frequencies, only the time-independent part of 
this sum is relevant in the above equations, leading to: 



id t <p 2 (t) ~ — 
ft 



(14) 
(15) 



where J e g is given by Eq. ([3]). 

The effective tunnelling can be both zero and negative. 
Note that sometimes the physics is described by a zero [5] 
whereas in other cases the dynamics is governed by the 
other Bessel functions |15j . Equations (14) and (15) yield 
the effective, time-independent Hamiltonian: 



eff 



(4, 



! 2 ^ ^2^1 



(16) 



Repeating the above derivation for the case of an optical 
lattice (see Appendix) gives the effective Hamiltonian ^ 
for non-interacting particles. It remains valid even for 
interacting particles, as long as the interaction is not too 
large. 

Strictly speaking, in the periodically shaken case there 
is no lowest energy: energies are only defined up to in- 
teger multiples of hu. Thus, the quasi-energies for high 
frequencies read: 



€ = ± J eff + jflW , 



(17) 



with integer j. 

For negative effective tunnelling, the ground state wave 
function in the effective model is antisymmetric, which 
would be impossible without shaking. Despite the infinite 



number of quasi-energies, there are only two Floquet- 
states for the case discussed here: 



± ^« t / h ( exp[^ sin(w< + a)] 



^2 V ± ex PH£ sin M + «)] I (18) 

In order to calculate the time-averaged energy in the 
time-dependent Hamiltonian, 



m))± = \ 



(19) 



one needs 



(i)±(t)\fi\i>±(t)) = -hfj,co8(ut + a)4>*(t)il>i(t) 
+h/icos(u)t + ajip^i^^it) 

-.mitWiit) ■ (20) 

The second step involves the time-average over one os- 
cillation period, which removes the parts proportional to 
ip*(t)ipj(t), j = 1,2 and leaves: 



(#(*)>- 



1 

T 



T 



= T^eff , 



(21) 



which uses Eq. ( 13 ). Note that the only approximation in 



the above equation involves the high-frequency approxi- 
mation for the Floquet-states ip±(t). 



IV. SUDDEN CHANGES OF THE 
HAMILTONIAN 

A. Forcing time-reversal via changes to the shaking 
amplitude 

If the Hamiltonian is changed quasi-instantaneously, 
this may involve changing, for example, either the ampli- 
tude or the phase or both at t = to: 



a 



Mi 
• a. 



(22) 
(23) 



To understand the effect of sudden changes of the Hamil- 
tonian, the example of a single particle in a double well 
is again a good starting point. In this case, the eigen- 
functions will be given by 
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e ±iJ cff t/R 


/ r 

( exp 1 




V ± ex P [ 


e ±i.7 off t/R | 


exp i 


V2 | 


± exp 



\{uit + a)] 



-i^sin(a;t + a) 



t < t 



t > t Q 



(24) 



Choosing the example £j = § shows that a particle which 
is initially in the V>+-state might be in the tp_ -state at 
t > to if the jump takes place for sm(ujto+a) = 1 and if a 
is chosen such that sin(wto + 5) = 0. Thus, although the 
wave function is continuous at t — to, the average energy 
can change even if the effective Hamiltonian is the same 
before and after the jump (cf. [27l 139] ). 

In the following, a quasi-instantaneous change of the 
Hamiltonian is used to switch the sign of the time- 
independent Hamiltonian ([2]). The sign can be switched 
by quasi-instantaneously changing both the tunnelling 
term by switching the shaking amplitude, e.g. pQ, 



Jo(1.692) ~ 0.403 and 
Jo(3.832) ~ -0.403 , 



(25) 
(26) 



and the sign of the interaction via a Feshbach- 
resonance l40l. 



#ideal = 



-H eS : < t < T 
-H oS : t > t 



U 



(27) 



The corresponding unitary time-evolution is given by 
C/(0,t) = exp[i(r — 2t q )H cS / J] for r > r with perfect 
return to the initial state at r = 2to- The turning point 
To however has to be chosen with care: only by taking 
To close to the maximum of the shaking, can unwanted 
excitations be excluded (cf., e.g., Ref. [27]). I n the fol- 
lowing, the time-reversal is demonstrated by numerically 
solving the full, time-dependent Hamiltonian ([I]) using 
the Shampine-Gordon routine |41j . 



B. Complex J c s induced via phase-jumps 



As the ansatz in Eqs. ^ and (10 1 is not unique, one 
can also derive effective equations for which J c g becomes 
complex [T51 [501 (35- Choosing, e.g., complex phase- 
factors such that the exponential in the ansatz is 1, i.e., 



ipi{t) = <fii(t) exp i— sa\(uot + a) — i— sin(a) 
L uj uj 

ip2{t) — f2{t) exp — i— sm(uit + a) + i— sin(a) 



would replace the effective Hamiltonian ( 16 ) by 



Gefi = — Ji 



cir 



— I— sin 



(Q) 4cO (28) 



where J c g still is the real number defined in Eq. ([3]). 



The corresponding 2 x 2-matrix has the same eigen- 
values but different eigenvectors than Eq. (16). How- 



ever, the full time-dependent solution is again identical 
to Eq. (18). Thus, as before, the ground state does not 
depend on the phase of the driving. Choosing 



i — sin(a) 



1 



(29) 



rather than 1 would, however, give the impression that 
the sign of the effective Hamiltonian has changed and 
ground and excited s tates h ad swapped their places (this 
was the aim of Sec. IV A). Using the full Hamiltonian 
shows that is not the case. In the following, the jump 
is modelled at t = under the assumption that before 
the jump one has a = and J e g > 0. Introducing 
the phase-jump (29) quasi-instantaneously changes the 
Floquet-states 



|V>-(f)>a=0 



\lp-(t))a?0 
|V>+(*))a#0 



(30) 
(31) 



where a ^ refers to the value defined by Eq. (29) but 



leaves the time-averaged energies ( 19 1 unchanged 



(H(t)) + = const. 
(H(t))_ = const. 



(32) 
(33) 



Contrary to this, the sudden change of the amplitude dis- 



cussed in Sec. IV A leaves the wave functions unchanged 
but quasi-instantaneously changes the sign of the average 
energy. This leads, in particular, to having a ground state 
which is antisymmetric. Another example of a related 
system with a counter-intuitive ground state is discussed 
in [32]. As the ground-state properties do not change 
by the jump (29), this paper avoids calling Eq. (28) an 



effective Hamiltonian. 

As long as one is primarily interested in the dynamics 
(rather than, e.g., the ground state) it would neverthe- 
less be perfectly valid to use such equations; the decision 
if phase-jumps or jumps in the shaking amplitude are 
preferable is a question of available technical laboratory 
resources and thus can not be discussed by the present 
paper. While this paper focuses on sudden changes for 
which sin(a) = 0, discussing sin(a) ^ can be easier for 
Eq. p8| rather than Eq. ( 16 ) as the latter will involve 



complicated changes in the effective wave function at the 
time of the jump. 
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(a) 




-14 



14 




t/jt 



FIG. 3. (Colour online) Time-reversal of the quantum dynam- 
ics for two interacting particles in a periodically shaken opti- 
cal lattice, (a) Two-dimensional projection of the probability 
density (to the power 1/4, in arbitrary units) if both parti- 
cles initially are close to the middle of the lattice (U — 1.6J, 
fiw — 20J, 2/i/oj = 1.692 J), (b) After applying coherently 
controlled time-reversal at r = to, the "echo" is visible near 
t ss 67r (to = 3n; all parameters as in a except for r > to: 
2fi/u = 3.832 J and NU = -1.6J). 



V. RESULTS 

Figure [3] shows the spreading of the probability den- 
sity for two interacting particles initially localised near 
the centre of a periodically shaken optical lattice. The 
initial state was chosen to be the ground state of a three- 
lattice-site version of the effective Hamiltonian ^ where 
the potential energy of the two outer wells was increased 
by 10 J compared to the middle well, thus mimicking the 
strong harmonic confinement of a Bose-Einstein conden- 
sate released into such a shaken lattice in Ref. [13]. To 
visualise the probability density (which is a single num- 
ber for each lattice site), a Gaussian density profile of the 
Wannier-function was modelled. 

The dynamics displayed in Fig. [3] (a) are similar to the 
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FIG. 4. (Colour online) Time-reversal in a per iodically shaken 
optical lattice (a) The return probability ( |34[ ) at r = 2ro as a 
function of interaction. It is only large if both interaction and 
amplitude are changed (red/triangles). Switching only the 
shaking amplitude (blue/circles) as expected only works for 
non-interacting particles, switching only the interaction (ma- 
genta/crosses) may work for high interactions, (b) The return 
probability as a function of switching time has a broader dis- 
tribution for low frequencies (blue/crosses, ftu = 10 J) than 
for high frequencies (red/circles, htu — 20J). All other pa- 
rameters are as in panel (b). (c) The return probability as a 
function of shaking frequency decreases for lower angular fre- 
quencies (all other para met ers as in b ). While in the high fre- 
quency regime [cf. Eqs. ( 18 1 and ( A9 1] instantaneously chang- 
ing the Hamiltonian should be equally reliable for any zero of 
sin(arf+a); for lower frequencies multiples of 2-zr lead to better 
return probabilities than ir, 3n, . . . (cf. [27]). 



text-book case of a free particle in one-dimension [33J . At 
first glance, the main difference to the classical diffusion 
equation seems to be the rate at which the probability 
density expands (cf. [34]). The time-reversal visible in 
Fig. [3] (b) demonstrates deeper differences between clas- 
sical and quantum physics. To quantify the quality of 
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time-reversal, the "return probability", i.e., the proba- 
bility to return to the original state, 



Appendix A: Effective Hamiltonian for a single 
particle in a tight-binding lattice 



P=|(*(t = 0)|*(t = 2to))| 2 , 



(34) 



is used (Fig. [4]). This figure shows that, as for the case 
of the Bose-Einstein condensate in a periodically shaken 
double well pQ, the time-reversal is not too sensitive to 
the precise choice of experimentally relevant parameters 
like interaction [Fig. [4] (a)], point of switching [Fig. 0(b)] 
and shaking frequency [Fig. [4] (c)] . 

The high probability to return to the initial state ob- 
served in the numerics offers an ideal method to dis- 
tinguish between statistical mixtures and pure quantum 
states (cf. pQ). For larger particle numbers, experiments 
could observe the "echo" by studying the width of the 
wave function with the existing experimental setups of 
RefsHDH31. 



VI. CONCLUSION 

To conclude, effective time-reversal was numerically in- 
duced via periodic shaking to an optical lattice for two 
indistinguishable interacting bosons. The time-reversal 
was ideally described by an effective Hamiltonian and 
subsequently tested numerically by using the full, time- 
dependent model. 

Time-reversal was realised by switching the sign of the 
effective Hamiltonian by quasi-instantaneously changing 
the strength of the shaking. At the same time, the 
scheme suggests a change of the sign of the interaction 
via a Feshbach resonance. The Floquet-states remain un- 
changed by this; however their time-averaged mean en- 
ergy switches its sign. This reflects the change of the sign 
of the effective Hamiltonian and supports the claim that 
the effective equations used to describe the dynamics can 
be labelled as an effective Hamiltonian. 

The high probability to return to the initial state shows 
that the time-reversal can be used to distinguish quan- 
tum mechanics from statistical mixtures induced by de- 
coherence. 

The numerical simulations indicate that the effective 
time-reversal should be realisable with present day tech- 
nology. To further optimise suitable time-reversal in an 
experiment, the feed-back loops of optimal control theory 
could be used EH SSI. 
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The Schrodinger equation for a single particle in a 
shaken lattice reads: 



idi^e(t) = 2i^cos{uit + a)ipe(t) 

- ^/+i(t) + y*-i(t)] 

with integer I. Choosing the ansatz 
where 



(Al) 



A(t) = exp 



2 M • , 
-i — sm(wt + a) 



(A2) 



(A3) 



leads to 
idnpe(t) 



J 



[A(i) mi (t) + [A(t)]-V/-i(*)]<A4) 



Within the high-frequency regime, one gets: 

^cff 



idnpi{t) 



(A5) 



where J e g is again given by Eq. (§. 

Thus, the effective Hamiltonian now reads: 



oo 



(A6) 



The eigenfunctions of this effective Hamiltonian are 
known: 

\<p k (t)) = exp (-i^-t) £ exp(ikdt)\e) , (A7) 



where d is the lattice constant, k the quasi-momentum 
and the energy is 



e(k) = — 2 J c ff cos(kd) 



(A8) 



The notation e(k) was chosen to indicate that it will be- 
come the quasi-energy of the time-dependent problem. 

Within the high-frequency approximation, the 
Floquet-states are now: 

\Mt)) = exp (-i^-t) [A(t)] e eMikd£)\e). 

(A9) 

In order to calculate the mean energy, time-averaged over 
one oscillation period T, where ujT = 2tt (i.e., the energy 
associated with this Floquet-state) we can replace the 
sum X)fc-oo by SiLi an d take the limit N — > oo 
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at the end. Thus, with this normalised version of | V'fc (*) ) Without any further approximation this leads to: 
defined on N lattice sites we can write: 

E k = ^£dt(Mt)\mMt)) (MO) E- k = -2JJ (^)cos(kd) (All) 

i N { i r T = — ^ eff cos (^^) 

r 1 " 3 ^ Jj \ t / dt 2£hfl C0S ^ + ^ = £ ( fc ) • ( Al2 ) 

o i I 



TV 

£=1 



T 

dt J exp(ikd) + [A(t) expiikd)}^ 1 ) 
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